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Abstract 


The radar return powers from a three-frequency radar, with center frequency at 22.235 
GHz and upper and lower frequencies chosen with equal water vapor absorption 
coefficients, can be used to estimate water vapor density and parameters of the 
precipitation. A linear combination of differential measurements between the center and 
lower frequencies on one hand and the upper and lower frequencies on the other provide 
an estimate of differential water vapor absorption. The coupling between the 
precipitation and water vapor estimates is generally weak but increases with bandwidth 
and the amount of non-Rayleigh scattering of the hydrometeors. The coupling leads to 
biases in the estimates of water vapor absorption that are related primarily to the phase 
state and the median mass diameter of the hydrometeors. For a down-looking radar, 
path-averaged estimates of water vapor absorption are possible under rain-free as well as 
raining conditions by using the surface returns at the three frequencies. Simulations of 
the water vapor attenuation retrieval show that the largest source of error typically arises 
from the variance in the measured radar return powers. Although the error can be 
mitigated by a combination of a high pulse repetition frequency, pulse compression, and 
averaging in range and time, the radar receiver must be stable over the averaging period. 
For fractional bandwidths of 20% or less, the potential exists for simultaneous 
measurements at the three frequencies with a single antenna and transceiver, thereby 
significantly reducing the cost and mass of the system. 
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1. Introduction 


Many remote sensing techniques are being used experimentally or operationally to 
estimate atmospheric water vapor. These include radiometers at microwave [Westwater 
and Han, 2001; Solheim et al., 1998; Schultz et al., 1993; Schlussel and Bauer, 1993], 

millimeter wave [Rosenkranz, 2001; Wang et al., 1995; Wilheit, 1990] and IR [Schmetz 

\ 

and Turpeinen, 1998; Hagen et al., 2004]; Raman and differential absorption (DIAL) 
lidars [Whiteman, 2003; Behrendt et al., 2002; Brassington, 1982]; and GPS-based 
techniques [Bevis et al., 1992; Alber et al., 1997]. With a few possible exceptions [e.g., 
Liljegren, 2004] most of the instruments and retrieval methods are not applicable in the 
presence of rain. 

Although little work seems to have been done on water vapor estimation using radar, an 
exception is the work of Tian et al. (2004) who have analyzed dual-frequency (10 and 94 
GHz) airborne Doppler radar data. By deriving the hydrometeor size distribution from 
the Doppler velocities and then modifying the radar reflectivity factors, Z, to account for 
Mie scattering and attenuation effects, the attenuation from cloud and gases can be 
inferred from the difference in the modified Z values at the two frequencies. In this paper, 
a three-frequency radar is studied where one of the frequencies is taken at the 22.235 
GHz line center with the others chosen at a lower and higher frequency with equal water- 
vapor absorption coefficients. Fractional bandwidths, defined as the ratio of the 
difference between the upper and lower frequencies to the center frequency, are 
considered at 0.1, 0.2 and 0.3. For an up-looking, ground-based geometry, retrievals of 
water-vapor absorption are possible only with precipitation present. For a down-looking 
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aircraft or spacecraft, path-integrated retrievals are possible in both rain and rain-free 
regions by using the surface as reference target. In the presence of precipitation, 
equations can be derived as functions of the radar return powers to any range at which 
backscattered powers are detected. The equations for precipitation and water-vapor 
attenuation are approximately separable in the sense that for small bandwidths, the 
equations for rain parameters and water-vapor absorption are weakly coupled. 

An error analysis of the approach shows that the biases in the method are relatively small 
in the rain and snow for percent bandwidths less than about 20% but become significant 
at bandwidths of 30%. These biases are related to non-Rayleigh scattering effects of the 
hydrometeors and are determined primarily by phase state and median mass diameter. 
The bias in the estimate of water vapor density is strongest in the layer of mixed phase 
hydrometeors associated with the melting layer. The largest source of variability in the 
estimates arises from signal fluctuations that can easily be as large as the differential 
signal itself. While the requirement for large numbers of independent samples (to reduce 
the signal variance) places demands on the radar design and operation, the relatively 
small frequency separation may allow a ‘differential-frequency’ implementation where a 
radar using a single antenna and transceiver can measure the three-frequency returns 
simultaneously. 

Equations for estimates of water vapor absorption and precipitation attenuation are given 
in section 2, followed by a simulation and error analysis of the water vapor absorption 
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and density estimates in section 3. Discussions on the algorithm and radar 
implementation are given in section 4. 

2. Equations for precipitation and water-vapor path-attenuation 

To simplify the final form of the equations, the following conventions are used. The 
‘measured’ radar reflectivity factor, Z m , at radar frequency f and range r is defined in 
terms of the radar return power, P r , by: 

Z m (f,r) = P r {f,r)r 2 l{c\Kj\) (1) 

where C is the radar constant and |K W | 2 is the dielectric factor of water which, by 
convention, is taken to be equal to its approximate value (0.93) for frequencies between 3 
GHz and 10 GHz and for temperatures between 0 C and 20 C [Battan, 1973]. The 
unattenuated radar reflectivity factor, or simply radar reflectivity factor, Z, is related to 
by: 

r 

Z(f, r) = Z m (/, r) exp{c \[k p (/, 5 ) + k c (f,s) + k v (/, *)]<&} (2) 

0 

where k p , kc, and k v are the specific attenuations from precipitation, cloud water, and 
water vapor respectively, and where the precipitation may include rain, snow and mixed- 
phase hydrometeors. More generally, contributions from oxygen and cloud ice should be 
added to (2). If the units of k are taken to be dB km' 1 then c=0.2xlnl0; if the units of k 
are in km' 1 then c=2. In this paper, we use the former. It is also convenient to define 2- 
way differential attenuations (dB) to range r for the precipitation arid cloud, A pc , and for 
the water vapor, A v , by: 
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( 3 ) 


A pc {fj,f-,r) = 2 \{[k p (fj,s) + *„(/,, 4)] - [*,(/„ 5) + * e (/„4)]}& 

0 

r 

AifpUr) = 2 j[* v (/,,s) - W/.*)]* (4) 

0 


We define 10 logio of Z and Z m by: 
£„(./>)= I Olog 10 Z m {f,r) 

Z(/,r)=101og 10 Z(/,r) 


(5) 


Throughout the paper we assume that measurements of P r , or equivalently Z m , are made 
at frequencies (fi, f, f u ) where f is taken at the center of the water-vapor absorption line 
at 22.235 GHz and (fj, f u ) are chosen such that fi < f c < f u and 


Kif,) = Kifu) (6) 

Using the water vapor model described by Ulaby et al. [1981] and Waters [1976], the 
values of (fi, f u ) that satisfy (6) are graphed in Fig. 1 as a function of the fractional 
bandwidth, Af, which we define by: 


A f = (fu~fl)lfc 


(7) 


If Af is specified, then (fi, f u ) can be found from the following approximations: 
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( 8 ) 


// =/ e -10.936A/ + 5.115A/ 2 
7>/ e+ 11.3A/+5.114A/ 2 (9) 

As a consequence of the skewed shape of the absorption line, k v (f c +5f) > k v (f c -5f) for 5f 
>0. For the calculations given later in the paper, percent bandwidths of 10%, 20% and 
30% are used where the frequency pairs, in GHz, are given, respectively, by: (f, f u ) = 
(21.248, 23.420), (20.246, 24.694), (19.409, 26.079). It is also worth mentioning that for 
percent bandwidths of 20% or less a “differential-frequency” implementation of the radar 
might be feasible. By this term, we mean that a single radar antenna and transceiver 
would provide measurements of the radar returns at the three frequencies . We will return 
to this issue later in the paper. 

Taking the difference of measured reflectivity factors (in dB) at the upper and lower 
frequencies, Z m (f u ,r)-Z m (f,,r) , and using (6) and the above definitions gives 

Z m (fu’ r )-Z m (f l ,r) = Z{f u ,r)-Z{f l , r )- A pc {f u ,f[\r) (10) 

Similarly, the difference between the measured reflectivity factors (in dB) at the center 
and lower frequencies, Z m (f c ,r)-Z m (f t ,r) , can be written: 


Z m <J e ,r)-Z m (fi,r) = Z(f c , r) - Z(/ ; , r) - (A pc (f c J,;r) + A v (f c , f t ; r)) 


(ID 
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To obtain an estimate of the 2-way differential water vapor absorption, 4X/o//) » that is 
approximately independent of the precipitation and cloud attenuation, we assume that 

A pc (f c ,fi\r) = yA pc (J u ,f l \r) (12) 

Eq. (12) states that the differential attenuation arising from precipitation and cloud 
between frequencies (f c , fi) is equal to some fraction y of the differential attenuation 
between frequencies (f u , fi). If the precipitation and cloud attenuation is directly 
proportional to frequency, y is simply the ratio of frequency differences: 

/ R ay=(fc -m/u-f,) (13) 

and (12) is satisfied exactly. However, because of the frequency dependence of the 
dielectric constant of water and non-Rayleigh scattering, y depends on the hydrometeor 
size distribution so that (12) will generally be incorrect. As will be shown in the results 
below, the magnitude of the error depends on the nature of the backscattering medium 
(rain, snow or mixed phase), the parameters of the hydrometeor size distribution, the 
fractional bandwidth of the radar frequencies that are used, and the choice of y. 

Multiplying (10) by y and subtracting (1 1) from the resulting equation yields: 

r Z m (./»+( l-r)Z m (f l ,r)-ZMc,r) = E l +E 2 +A v (J c J l ;r) (14) 

where 
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El = rZ(f u ,r) + (1 - r)Z(f t ,r) - Z(f c ,r) 


(15) 


E 2 =A pc (f c ,f l ;r)- r A pc (J u ,f l ;r ) (16) 

The bias term E 2 was discussed above where it was noted that E 2 =0 if (12) holds. The 
bias term Ei is a function of the reflectivity factors at the 3 frequencies. Unlike the 
measured reflectivity factors on the left-hand side of (14), they are not measurable and 
not easily estimated unless the precipitation retrieval problem is solved in parallel. 

Similar considerations apply to this term as apply to (12): if the backscattering is 
Rayleigh, the reflectivity factors are approximately frequency independent and Ej £ 0. 

To estimate A,(f c ,f,) from measured data, we assume that the error terms are zero so that 

A y ( f c ,f/',r) = yZ m (/„ , r) + (1 - y)Z m (f, ,r ) - Z m (f c ,r) (17) 

Before presenting an error analysis of (17), we derive equations similar to (10) and (14) 
by using the backscattered powers from the surface. The equations below are obviously 
relevant only to a down-looking geometry. On the other hand, unlike (10) and (14), they 
are applicable in rain-free as well as raining conditions. The measured or apparent radar 
return power from the surface, P sm , can be written as [Kozu, 1995]: 

rs 

p sm (f A s ) = [C, (/)o-° (/) / gfo)] exp {-0.2 In 1 0 j[£ p (f,s) + k c (/, s) + k v (f, s)] ds} ( 1 8) 

0 

where a 0 is the normalized radar cross section of the surface (unitless) and C s is the radar 
constant for surface scattering which depends on incidence angle, antenna gain, and 
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frequency. The function g(r s ) represents the range dependence which is equal tor/ for 
beam-limited (nadir and near-nadir incidence) and r/ for pulse-limited conditions (off- 
nadir incidence). The radar return power from the surface, P s , which would be measured 
in the absence of atmospheric attenuation, can be written: 

P s (f,r s ) = C s (f,r s )<j\f)lg(r s ) (19) 

As before, we define the following quantities: 

P s (/) = 10 log 10 P s (f) 

Psm(f) = 10 log, o P !m (f) (20) 

a°(/) = 101og ]0 CT 0 (/) 


Comparing (18)-(20) with (1), (2) and (5) shows that for the case of surface scattering, 
(10) and (14) become, respectively: 


Psm(fu’ r s)~ Psm(fl> r s) = Ps(f«> r s)-Ps(fb r s)- A pc(fu’fl'’ r s) 


rPsn,(fu’ r s) + Q- r)Psm(fl’ r s)- Psm(fc’ r s) = E ls + E 2 s + MfcJb r s) 


( 21 ) 

( 22 ) 


where 


E u = rPs ifu ,r s ) + (\-y)P s (/, ,r s )-P s (f c ,r s ) 

E 2s= E 2( r = r s) = A pc ( f c r A pc ( fu Jb r s) 


(23) 

(24) 
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One difference between the rain and surface scattering equations above is that the Ei s 
error depends on the nature of backscattering not from the hydrometeors, as in Ei, but 
from the surface. This becomes apparent if (23) is written as: 

E u = [y°\f u ,r s ) + {\-y)a\f l ,r s )-5\f c ,r s ))+... (25) 

where the notation is used to indicate additional terms that are functions of the 
radar constants. Since these constants are known, or presumed to be known, they can be 
accounted for in an offset term. A comparison of the bracketed term in (25) with (15) 
shows that the approximation for the surface scattering is analogous to that used for the 
precipitation, with the normalized surface cross sections replacing the radar reflectivity 
factors. If the surface cross section is constant over the frequency span from fi to f u the 
bracketed term in (25) is zero for any y. Moreover, if y is set to yRay, given by (13), then it 
can be shown that this term is zero if <?° changes linearly with frequency. Another 
difference between Ei and Ei s is that an estimate for the latter quantity can be obtained 
from surface scattering measurements in regions of low path-integrated water vapor. 

This is analogous to the ‘surface-reference technique’ where measurements of the surface 
return in rain-free areas are used to estimate path-attenuation in the presence of rain [Tian 
et al., 2002; Meneghini et al., 2004]. In particular, an estimate for the differential path- 
attenuation from cloud and hydrometeors follows directly from (21) where the term 
? s (f u ,r s )-P s (fi,r s ) is approximated by surface returns measured in clear conditions. 

Since the main subject of the paper is the feasibility of height-profiled water vapor 
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estimates in rain, we focus the error analysis on the range-dependent estimate of 
differential water vapor attenuation given by (17). 

3. Error analysis of water-vapor profiling retrievals 
3a. Description of the simulation 

To analyze errors in the estimate of differential water vapor given by (17) we construct a 
simple stratiform storm model derived from disdrometer-measured raindrop size 
distributions [Meneghini et al., 2003]. For each size distribution we compute the median 
mass diameter, Do (mm) and number concentration, N t (m' 3 ). To complete the 
specification, we assume that the drop diameter distribution, N(D) (mm 1 m' 3 ), follows a 
gamma distribution with a fixed shape parameter, p, equal to 2 [Ulbrich, 1983]: 

N(D) = N 0 D M exp[-(3.67 + ju)D/D 0 \ (26) 

where 

N 0 = N, (3.67 + M y +1 /(Df'np + 1)) (27) 

and where T is the gamma function given by: 

oo 

T(x)= jt xA e-'dt (28) 

o 

The storm structure is assumed to be composed of a layer of snow, from 4 to 5 km above 
the surface, a transition region of mixed-phase hydrometeors, from about 3.5 to 4 km and 
a rain layer from about 3.5 km to the surface. The fractional melt water in the melting 
layer is prescribed as a function of the initial snow size, mass density and distance below 
the 0° isotherm [Yokoyama and Tanaka, 1984]. The effective dielectric constant of the 
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mixed phase particles is given by the effective medium approximation [Bohren and 
Battan, 1980]. For the results shown here, the snow mass density is fixed and equal to 
0.2 g cm' 3 . The location and density of the cloud water are also included although this 
contribution is not critical to the behavior of the water vapor retrievals. The temperature 
lapse rate is taken to be constant and equal to 6 C/km so that with the 0° isotherm at 4 km, 
the surface temperature is 24 C. To focus on the bias and random errors of the estimate 
in this set of examples, the relative humidity is taken to be constant from the surface to a 
height of 5 km and equal to 80%. From this simple storm model, approximately 400 
range profiles of Z m at the three frequencies can be generated, where each profile is 
derived from measured raindrop size distribution parameters and consists of the returns 
from the snow, mixed phase and liquid hydrometeors, and includes effects of absorption 
from cloud water, water vapor and molecular oxygen. Note also that the range resolution 
is taken to be 125 m so that the profile consists of 40 gates or bins. The simulated radar 
measurements at f = fi, fi, f u are computed from the following set of equations: 


Z m (J*r) = c z [ Jo- 6 (/, D)N(D, 5>/£>]exp[-0.2 In 1 0 J(*„ (/, s) +k c (f,s) + k v (f,s))ds] (29) 

0 0 

where the specific attenuations (dB km' 1 ) are given by: 

oo 

*,(/> 4.343 xl0~ 3 \a e (f,D)N(D,s)dD (31) 

0 

^(/) = [(4.343xlO- 2 x6#)/c 0 ]Im(-^ w )M c (32) 

& v (/) = 2 / 2 /? v (300/r) 3/2 yi((300/r)e _644/r /[(494.4-/ 2 ) 2 + 4/ 2 y 1 2 ] + 1.2xl0 -6 } (33) 

with 
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Cz =c 0 4 /(/V|KJ 2 ) 


(34) 


y x = 2.85(/7I01 3)(300 / T 1 ) 0 626 [1 + 0.0 1 8 / 7 / / />] 

where P is the pressure (millibar), T is temperature (K), p v is the water vapor density (g 
m' ), Co is the speed of light (mm s' ), m is the complex index of refraction of the (cloud) 
water, M c is the cloud water content (g m' 3 ), and <Tb(f, D), a e (f, D) are, respectively, the 
backscattering and extinction cross sections (mm 2 ) of a sphere of diameter D at frequency 
f. The expression for k v above is taken from results of Waters (1976) and Ulaby et al. 
(1981) for f < 100 GHz. It is important to note that in (33) the frequency is to be 
specified in GHz; in all other equations the frequency is to be expressed in Hz. 

Following Ulaby et al. (1981), the contribution from molecular oxygen is also included. 

The error analysis consists of computing the right-hand side of (17) using (29) - (36) and 
comparing the estimated 2-way differential vapor absorption profile with the assumed 
profile as calculated from (4), (33) and (36). To compute the right-hand side of (17), a 
value of y is required. Although the approximation from (13) can be used, it was found 
by trial and error that somewhat smaller values reduce the bias. For the results shown in 
the paper, we use y=0.44 (10% bandwidth), y=0.42 (20%), and y=0.39 (30%) instead of 
the results from (13) which give: yRa y =0.455 (10% bandwidth), yRa y =0.447 (20%), and 
YRay =0.424 (30%). An exception to these assumptions is the case presented in Fig. 2, 
explained below, where we have chosen y=0.39 for the 20% bandwidth case to emphasize 


(35) 

(36) 
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the bias in the estimate; i.e., the bias is qualitatively the same but smaller when the value 
y=0.42 is used. 

Apart from the bias errors, Ei and E 2 , described in the previous section, an important 
additional source of error arises from the finite number of independent samples used to 
estimate the radar return power and Z m . For a square-law detector, where the output of 
each sample is proportional to the return power, the deterministic quantity Z m is replaced 
by ( Z m + g ) , where, for the number of samples, n, greater than about 1 0, the random 

variable g can be approximated by a zero mean Gaussian with variance [Bringi et al., 
1983]: 

n - 1 

var(g) = 4.343 2 {k 2 / 6 -^m ~ 2 } (37) 

m = 1 

For figures 2 through 8, n is taken to be a large number (64,000) so that the bias terms 
can be examined without large background variability from finite sampling. 

3b. Error estimate of differential water vapor attenuation 
The top panel of Fig. 2 shows simulations of Z m (f c ) over approximately 400 
‘observations’ (x-axis) over the 5 km storm height (y-axis). The results clearly show a 
bright-band, corresponding to returns from the partially melted snow, in the ranges just 
below the 0° isotherm at 4 km. Since the radar is assumed to be above the storm viewing 
along nadir, effects of attenuation are evident by the reduction in the measured 
reflectivity factor as the penetration depth increases. Results for the measured 
differential reflectivity for the upper and lower frequencies, Z m (/„)-Z m (/ ; ), are shown 
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in the second panel from the top. A 20% bandwidth has been chosen so that 
f u = 24.694 GHz and f, = 20.24 6GHz . The magnitude of Z m (/„ ) - Z m (/, ) depends on the 
drop size distribution, phase of the hydrometeors (solid, liquid or partially melted), and 
the cumulative attenuation by precipitation and cloud. However, it is independent of 
water vapor since k v (f u ) = k v (f,). 

The bottom two panels of Fig. 2 show the estimated and assumed values of the 2-way 
differential water vapor absorption, A t (f c ,f l ;r) , where the estimated value is computed 
from the right-hand side of (17). Despite the large number of independent samples 
(64,000), variations caused by the finite number of samples are still evident. 
Superimposed on the fluctuations are biases that are seen to be well correlated with 
decreases in the Z m (/„)-Z m (/,) field shown in the panel above; as already noted, the 
bias terms have been exaggerated by choosing y=0.39 for this plot. 

To understand how the bias in A v {f c ,f l \r) depends on the parameters of the precipitation, 
it is useful to examine some of the relevant quantities as functions of the median mass 
diameter. Do. The results in Figs. 3 and 4 are taken from the lowest gate, just above the 
surface so that the attenuations are approximately the same as the total path attenuations. 
Shown in the top panel of Fig. 3 are values of Z m (f u )-Z m (f,) versus Do. Variations in 
this quantity originate from two sources: non-Rayleigh scattering and attenuation effects 
from the precipitation and cloud. The non-Rayleigh source of variability can be seen in 
the results for Z(/„)-Z(/,) shown in the middle panel of Fig. 3. Since the ‘shape’ 
parameter, p, of the size distribution has been fixed, Z(/„)-Z(/ / ) is a function only of Do 
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and independent of the number concentration. A scatter plot of the differential path 
attenuation from cloud and precipitation, A pc (f u ,f t \r) , is shown in the bottom panel. 

Note that this quantity depends on N t as well as Do as seen by the variability in this 
quantity for fixed Do- Since the cloud liquid water was taken to be zero for this 
simulation, the attenuation arises entirely from the hydrometeors. The addition of an 
integrated cloud water content of 1 kg/m 2 for all cases yields an increase in 
Apc(f u ,f,; r ) of about 0.4 dB (not shown) but has a negligible effect on the bias errors 

shown in the lower two panels of Fig. 4. 

The top panel of Fig. 4 shows the estimated (X data points) and true values (solid line) of 
the differential water vapor absorption, A v (f c ,f,;r) . The variability in the estimates for 
constant Do is the result of finite sampling and reduces to zero as n goes to infinity. In 
addition to the variability is a bias which is negative for Do values up to about 2 mm and 
slightly positive above 2 mm. This bias is the sum of Ei and E 2 displayed in the bottom 
two panels of Fig. 4. The behavior ofEi is determined almost exclusively by the median 
mass diameter; on the other hand, variability at a fixed Do, indicating dependence on 
number concentration, can be seen in the E 2 term for large D 0 values. The results imply 
that if the rain estimation problem can be solved in parallel (or iteratively) with the water 
absorption estimate, the Do (and N t ) values, as estimated from the precipitation 
algorithms can, in principle, be used to correct for biases in the water vapor retrieval. 
Although a quantitative description of algorithms of this type is beyond the scope of the 
paper, a qualitative description of a possible approach is outlined in section 4. 
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The results in Figs. 3 and 4 are taken from the range gate in rain at the bottom of the 
column. For the case shown in Figs. 5 and 6, the range is taken to be just above the 
melting layer at a snow depth of about 0.9 km. Comparison of the results of Fig. 5 and 
Fig. 3 show a number of obvious differences between the two cases. The differential 
reflectivity factor in snow ,Z(f u )-Z(f t ) , shown in the middle panel of Fig. 5, is a 
monotonically decreasing function of Do (the median mass diameter of the melted snow); 
moreover, as the attenuation through snow is small (bottom panel), the behavior of the 
measured differential reflectivity factor, Z m ( f u ) - Z m (/, ) , shown in the top panel, is 
primarily determined by Z(f u )-Z(f,) . 

The estimated and true values of the water vapor absorption to this range are shown in the 
top panel of Fig. 6. Notice that the magnitude of the fluctuations is essentially the same 
as before but because the water vapor absorption to this range is small, the relative error 
caused by the finite sampling is much higher. On the other hand, the biases in the 
estimate are small. This can be seen in the bottom panels where E 2 is approximately zero 
while Ei is small up to moderate Do values and positive for large Do. 

In the melting layer, the E 2 bias term remains relatively small; however, the Ei bias term 
is significant and changes from gate to gate. Plots of Ei versus Do are shown in Fig. 7 for 
four range gates within the melting layer. For gates 13 and 15, the sign of the bias is 
negative for small Do and positive for larger values. As will be shown later, this bias has 
a strong effect on the water vapor retrievals in and about the melting layer. 


18 



Before looking at statistics of the A v (f c ,f l ;r) estimates, it is instructive to view individual 
retrievals. In Fig. 8, five arbitrarily chosen profiles of A, (f c , f,\r) , from the results of Fig. 
1, are represented by the thin solid lines. The true (assumed) profile is represented by the 
heavy solid line. The results for the nearly 400 profiles are indicated by the + signs; that 
is, at each range gate, approximately 400 data points are displayed corresponding to the 
estimated values of A^{f c ,f l \r') at range r. The behavior of the individual profiles 
indicate that even for a very large number of independent samples, the fluctuations in 
A(fc>fi> r ) range can be substantial. Computing the differential specific absorption 
without smoothing in range or further averaging in space or time would clearly lead to 
large errors. 

The results shown in Fig. 8 are from the single simulation shown in Fig. 1. To obtain 
accurate estimates of the mean and standard deviation as a function of range, 200 
simulations of the retrievals were made. Note that from realization to realization, the set 
of drop size distributions is fixed and only the Z m are changed to simulate the effects of 
finite sampling. The results are shown in Figs. 9-11 for bandwidths of 20%, 30% and 
1 0%, respectively. Unlike the previous results, a more realistic number of independent 
samples (n=4000) is used. For each figure, images of the estimated and true values of the 
differential absorption, (f c ,f,;r) , for a single simulation are shown in the top panels; in 

the bottom panel are shown the true value (solid line) along with the mean (X) and twice 
the standard deviation of the estimate (vertical height of error bar) as a function of range. 
As in Fig. 8, gate 40 represents the range gate just above the surface. 
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The three sets of results show that the standard deviation is independent of range so that 
the relative error decreases with distance from the radar since the path absorption 
increases. In snow, the bias in f\r) is small. In the rain, however, the bias for the 
30% case (Fig. 10) is relatively large. The bias decreases in going to smaller bandwidths 
so that the 10% case (Fig. 1 1) shows the smallest bias. On the other hand, the fractional 
standard deviation (standard deviation divided by the mean) is largest for the 10% case 
and decreases in going to 20% and again to 30% bandwidth. 

3c. Error estimates of water vapor retrievals 

A range derivative of the differential path attenuation, A v {f c ,f l '), yields the differential 
specific water vapor absorption, k v (f c )-k v (f t ) . To reduce the variability in A(fc’fi) > a 
5-gate moving average is done, followed by a differencing over successive 5-gate 
intervals. After normalizing by the range difference, this procedure gives 
K(f c i r j)~k v (f l ,r j y,j = 3,..., 38 where estimates at the two lower-most and upper-most 

range gates are discarded. From the estimate of differential specific absorption and 
model temperature (with lapse rate as before of 6 C km’ 1 ) and pressure profiles, the water 
vapor density and relative humidity are computed. The ‘true’ temperature and pressure 
profiles are assumed to be Gaussian random variables with a mean given by the model 
profiles with standard deviations of 1 K and 2 millibar, respectively, as derived from the 
results of 200 soundings in the South China Sea. Note that the ‘true’ p v field is computed 
from the assumed RH and randomly varying temperature and pressure profiles. The 
simulated storm parameters are similar to those used previously with the exceptions that 
the integrated cloud water content is taken to be 1 kg/m 2 and the relative humidity is 
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assumed to be 100% above the 0° isotherm, decreasing linearly to 70% at the surface. 
Results are shown only for 20% bandwidth case with the number of samples, n, equal to 
1 6,000 samples. For n=64,000 the RMS errors in the water vapor and RH retrievals in 
the lowest 3 km are approximately halved while for n=4,000 the errors are approximately 
doubled. 

Shown in the top left panel of Fig. 13 are the estimated (dashed line) and assumed (solid 
line) mean values of A v (f c ,f l ) computed over the set of about 400 profiles. On the top 
right panel are the estimated (dashed line) and assumed mean values of 
K(fc)-Kifi) using the averaging procedures described above. Biases in f c ,f,) can be 
seen to increase just below the melting layer at 4 km and assume nearly a constant value 
from about 3 km down to the surface. Examination of the K(f c )-k v (f,) estimate reveals 
a large negative bias associated with scattering in the melting layer. This is directly 
related to the Ei bias term shown in Fig. 7. However, because of the need to smooth the 
data, the influence of this error extends into the snow above and the rain below. Standard 
deviations of 4,(f c >fi) and k v (f c )-k v (f,) are shown in the bottom panels. Since the 
assumed pressure and temperature profiles were randomized, the ‘true’ values of 

an d are also random with standard deviations given by the solid 

lines. The standard deviations of the estimated profiles of ^ (/<.,/,) and k v (f c )-k v (f,) are 
given by the dashed lines in the lower panels; the greater variability in estimated profiles 
is caused by the finite number of samples. As noted above, increasing n by a factor of 4 
reduces the estimated standard deviations by about a factor of 2. 
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Corresponding statistics for the water vapor density, p v , and RH are shown in Fig. 14 
where it can be seen that negative biases in k v (f c )-k v (f,) at and about the melting layer 
produce large negative biases in p v and RH over the same range. Although the vertical 
extent of the bias can be narrowed with an averaging interval less than five, this leads to 
higher standard deviations in the estimates. As in the previous plot, the standard 
deviations of the true and estimated fields are shown in the lower panels. 

The RMS errors for p v and RH are presented in the top panels of Fig. 15 and the 
normalized errors are shown in the bottom panels, i.e., the RMS error divided by the true 
value. For p v , the relative errors in the bottom 3 km are between 20% and 28%; above 3 
km, the errors become as large as 32%. Similar accuracies can be seen in RH, with 
relative errors between 20% and 28% below 3 km and larger errors above 3 km. If the 
assumed temperature and pressure are taken equal to the model values, the relative errors 
decrease by about 4%. Increasing the sample number to 64,000 decreases the relative 
errors in p v and RH to between 12% to 16% in the lowest 3km. In the melting layer, 
however, the maximum relative error is about 25%. Increasing the bandwidth to 30% 
decreases the standard deviations in k v (f c )-k v (f ,) » Pv, and RH in the lowest 3 km but 
also produces higher negative biases and relative errors in and about the melting layer. 
For the case of n=16,000 and 30% bandwidth, the standard deviations in p v and RH 
decrease to about 16% in the lowest 3 km but exhibit a maximum relative error of about 
36% in the melting layer. 

More sophisticated averaging and smoothing techniques might serve to reduce the RMS 
errors. As already pointed out, information on the median mass diameter of the 
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hydrometeors, particularly in the melting layer, would serve to reduce the bias. Further 
investigations, and, ultimately, an analysis of measured data, will be needed to address 
these issues. 

4. Discussion and Summary 

The approach used here has some similarities to the differential absorption lidar (DIAL) 
technique in that it uses frequencies on and offline center to estimate the strength of 
absorption. The fact that precipitation is the background scattering medium, however, 
implies that the differential attenuation by hydrometeors can easily be as large as water 
vapor absorption; moreover, non-Rayleigh backscattering effects at these frequencies can 
be comparable in magnitude to attenuation and absorption. By the use of 3-frequencies, 
we can take advantage of the fact that differential attenuation from precipitation and 
cloud is approximately an odd function with respect to the center frequency while the 
differential water vapor absorption is approximately an even function. 

In deriving an expression for the differential vapor absorption, the critical assumption is 
that the differential attenuation from cloud and precipitation between the upper and lower 
frequencies can be expressed as a fraction of the differential attenuation between the 
center and lower frequencies. Results of the simulations show that this assumption leads 
to biases in the estimate of water vapor density that become larger as the bandwidth 
increases. While the biases are small in snow, they are particularly strong in the melting 
layer and lead to large negative biases in the vapor density estimates in and around the 
melting layer. 
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Although the focus of the paper is estimation of water vapor, the approach also offers the 
potential of precipitation estimation. Because of the choice of lower and upper 
frequencies, the differential measured reflectivity factor, Z m (f u is function only 
of the characteristics of the precipitation and cloud and is independent of water vapor. If 
the differential path attenuation from cloud and precipitation can be estimated then 

follows directly from (10). But it is clear from the results of Fig. 3 and Fig. 
5 that from Z(/„)-Z(/ / )an estimate ofthe Do can be obtained. Moreover, an estimate of 
the number concentration N t can be obtained from the radar equation, (1). To start the 
procedure requires an initial or estimated path attenuation. One such estimate can be 
obtained from (21) by measuring the surface return powers in clear regions. However, as 
recently shown by Mardiana et al. (2004), the equations also can be solved iteratively 
without an independent path-attenuation estimate. In either case, the procedure yields 
estimates of the size distribution parameters along the profile. The Do and N t values can 
be used, in turn, to improve the estimate of the differential water vapor absorption by 
using the information to estimate the bias terms Ei and E 2 . A drawback to the procedure 
is that the equation for N t , as derived from the radar equation, is a function k v (fi) or k v (f u ). 
Although the term is usually small for the 20% and 30% bandwidth cases relative to the 
hydrometeor attenuation, it represents an additional error source in the precipitation 
retrieval problem. In principle, just as the hydrometeor size distribution parameters can 
be used to correct for biases in the water vapor retrieval, the water vapor retrieval can be 
used to account for this contribution in the precipitation retrieval. Whether iterating 
between solutions to the precipitation and water vapor equations will provide stable 
solutions is not clear, however. 
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Another way to view the approach is as a variation of the differential-frequency 
implementation. In this approach, a single antenna and transceiver are used to transmit 
and receive signals at more than one frequency. However, this requires a wide-band 
power amplifier and a wide-band antenna. Wide-band power amplifiers with bandwidths 
up to 20% are now available in some frequency bands [J. Carswell, personal 
communication]. For broad-band solid state amplifiers, with high duty cycles but low 
peak powers, pulse-compression techniques can be used to achieve fine range resolution. 
Averaging the data to a coarser vertical resolution may provide a sufficient number of 
independent samples to make the measurement technique feasible without excessive 
space or time averaging. Two other methods of increasing the effective number of 
independent samples are frequency agility and data “whitening” [Koivunen and Kostinski, 
1999; Torres and Zmic, 2003; Torres et al., 2004]. Another requirement of the radar 
would be well-matched beamwidths at the three frequencies. Since the total differential 
path absorption for a 20% bandwidth is on the order of 1 dB, any mismatches in the radar 
resolution volumes will have a strong effect on accuracy particularly in convective rain 
where vertical and horizontal gradients in the reflectivity field can be large. Hom-lens 
and parabolic antennas are inherently broad-band and should be capable of good 
performance over a 20% bandwidth. Nevertheless, detailed calculations would be needed 
to assess the degree of beam matching needed relative to the gradients in the reflectivity 
field. 

Acknowledgements : We wish to thank Professor Ramesh C. Srivastava for suggesting 
the use of multiple frequencies for the measurement of water vapor. 
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Figure Captions: 


Fig. 1 : Top: Specific absorption of water vapor (dB/km) for water vapor densities (g/m 3 ) 
of 8.75 (lowest), 13.15, 17.5 (highest) and the frequencies at which the absorption is 
equal for bandwidths of 10%, 20% and 30%. Bottom: Upper and lower frequencies 
giving equal water vapor absorptions versus percent bandwidth. Solid lines represent 
numerical values; dashed lines represent quadratic approximations given by (8) and (9). 

Fig. 2: Simulated fields of measured reflectivity factor at f c =22.235 GHz (top), 
differential measured reflectivity factor (second from top), estimated differential water 
vapor absorption (third from top) and assumed differential water vapor absorption 
(bottom). All units are in dB. (y=0.39). 

Fig. 3: Measured (top) and actual (center) differential reflectivity factor and differential 
path attenuation from precipitation (bottom) for a range gate near the surface versus the 
median drop diameter, Do. All units are in dB. 

Fig. 4: Top: Assumed (solid line) and estimated differential water vapor absorption (X) 
versus median raindrop diameter for a range gate near the surface; Center: Bias term Ei; 
Bottom: Bias term E 2 . y=0.42. 


32 


Fig. 5: Measured (top) and ‘actual’ (center) differential reflectivity factor and differential 
path attenuation (bottom) for a range gate in the snow just above the melting layer versus 
median drop diameter of melted particles, Do. 

Fig. 6: Top: Assumed (solid line) and estimated differential water vapor absorption (X) 
in snow versus median drop diameter of melted particles; Center: Bias term Ej; Bottom: 
Bias term E 2 . y=0.42. 

Fig. 7: Bias term Ei for four range gates within the melting layer (1 1, 13, 15, 17) versus 
median mass diameter of the equivolume raindrops. 

Fig. 8: Five profiles of estimated differential water vapor absorption (thin lines) and 
assumed differential absorption (thick line) versus range gate number (125 m) from storm 
top. Results from all estimated profiles are indicated by the ‘+” symbol. y=0.42. 

Fig. 9: Estimated (top) and assumed (center) differential water vapor absorption; bottom: 
assumed (solid line) and the mean (X) and twice the standard deviation (vertical bars) of 
the estimates versus range gate number. A bandwidth of 20% is assumed with y=0.42. 

Fig. 1 0: Same as Fig. 9 but for a bandwidth of 30%. y=0.39. 

Fig. 1 1 : Same as Fig. 10 but for a bandwidth of 10%. y=0.44. 
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Fig. 12. Top left: Mean profiles of estimated (dashed line) and assumed (solid) 
differential water vapor attenuation versus height. Top right: Mean profiles of estimated 
(dashed line) and assumed (solid) differential specific attenuation. Bottom left: Standard 
deviation of estimated (dashed line) and assumed (solid) differential water vapor 
attenuation. Bottom right: Standard deviation of estimated (dashed line) and assumed 
specific differential attenuation. n=T6,000 and 20% bandwidth. 

Fig. 13. Top left: Mean profiles of estimated (dashed line) and assumed (solid) water 
vapor density versus height. Top right: Mean profiles of estimated (dashed line) and 
assumed (solid) relative humidity, RH. Bottom left: Standard deviation of estimated 
(dashed line) and assumed (solid) water vapor density. Bottom right: Standard deviation 
of estimated (dashed line) relative humidity. n=l 6,000 and 20% bandwidth. 

Fig. 14. Top left: RMS error in estimate of water vapor density versus height. Top right: 
RMS error in estimate of relative humidity. Bottom left: Relative error in water vapor 
density estimate (RMS normalized by true value). Bottom right: Relative error in relative 
humidity estimate. n=l 6,000 and 20% bandwidth. 
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Fig.l. Top: Specific absorption of water vapor (dB/km) for water vapor densities (g/m ) 
of 8.75 (lowest), 13.15, 17.5 (highest) and the frequencies at which the absorption is 
equal for bandwidths of 10%, 20% and 30%. Bottom: Upper and lower frequencies 
giving equal water vapor absorptions versus percent bandwidth. Solid lines represent 
numerical values; dashed lines represent quadratic approximations given by (8) and (9). 
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Fig. 2. Simulated fields of measured reflectivity factor at f c =22.235 GHz (top), 
differential measured reflectivity factor (second from top), estimated differential water 
vapor absorption (third from top) and assumed differential water vapor absorption 
(bottom). All units are in dB. (y=0.39). 
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Fig.3. Measured (top) and actual (center) differential reflectivity factor and differential 
path attenuation from precipitation (bottom) for a range gate near the surface versus the 
median drop diameter, D 0 . All units are in dB. 
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Fig.4. Top: Assumed (solid line) and estimated differential water vapor absorption (X) 
versus median raindrop diameter for a range gate near the surface; Center: Bias term Ei; 
Bottom: Bias term E 2 versus median drop diameter. y=0.42. 
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Fig. 5. Measured (top), and ‘actual’ (center) differential reflectivity factor and differential 
path attenuation (bottom) for a range gate in the snow just above the melting layer versus 
median drop diameter of melted particles, Do. 


39 






0.0 


XX X»< X 


-o.i L 

0.5 


Fig.6. Top: Assumed (solid line 
snow versus median drop diam< 
Bias term E 2 . y=0.42. 
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Fig. 8. Five profiles of estimated differential water vapor absorption (thin lines) and 
assumed differential absorption (thick line) versus range gate number from storm top. 
The range gate spacing is 125 m. Results from all estimated profiles are indicated by the 
4 +” symbol. 
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Fig.9. Estimated (top) and assumed (center) differential water vapor absorption; bottom: 
assumed (solid line) and the mean (X) and twice the standard deviation (vertical bars) of 
the estimates versus range gate number. A bandwidth of 20% is assumed with y=0.42. 
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Fig.10. Same as Fig. 8 but for a bandwidth of 30%. y=0.39 
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Fig.l 1. Same as Fig. 8 but for a bandwidth of 10%. y=0.44. 
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Fig. 12. Top left: Mean profiles of estimated (dashed line) and assumed (solid) differential 
water vapor attenuation versus height. Top right: Mean profiles of estimated (dashed line) 
and assumed (solid) differential specific attenuation. Bottom left: Standard deviation of 
estimated (dashed line) and assumed (solid) differential water vapor attenuation. Bottom 
right: Standard deviation of estimated (dashed line) and assumed specific differential 
attenuation. n=l 6,000 and 20% bandwidth. 
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Fig. 13. Top left: Mean profiles of estimated (dashed line) and assumed (solid) water 
vapor density versus height. Top right: Mean profiles of estimated (dashed line) and 
assumed (solid) relative humidity, RH. Bottom left: Standard deviation of estimated 
(dashed line) and assumed (solid) water vapor density. Bottom right: Standard deviation 
of estimated (dashed line) relative humidity. n= 16,000 and 20% bandwidth. 
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Fig. 14. Top left: RMS error in estimate of water vapor density versus height. Top right: 
RMS error in estimate of relative humidity. Bottom left: Relative error in water vapor 
density estimate (RMS normalized by true value). Bottom right: Relative error in relative 
humidity estimate. n=l 6,000 and 20% bandwidth. 
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